function s=Area(A,B,C)
if length(A)==2
    AB=[B-A 0];
    BC=[C-B 0];
elseif length(A)==3
    AB=B-A;
    BC=C-B;
end
    Z=cross(AB,BC);
    s=1/2*norm(Z);
end